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Abstract 

We consider tlie real emission QCD correction to lieavy Higgs boson produc- 
tion via weak boson fusion in high energy pp cohisions. The 0{as) correc- 
tions are determined for the complete electroweak qq qqW~^W~ process. 
The presence of a third parton in the final state affects the formation of ra- 
pidity gaps only slightly. In particular, soft emission into the gap region is 
severely suppressed. Also, we investigate how the additional hard emission 
affects forward-jet-tagging and central-jet- vetoing efficiencies in the search for 
H W+W- l^vl'v decays. 
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I. INTRODUCTION 



One of the prime objectives of experiments at hadron supercolliders is the discovery of the 
Higgs boson and the subsequent investigation of its properties. If the Higgs boson is relatively 
heavy or if electroweak symmetry breaking is driven by some new strong interactions, then 
the study of Higgs boson production in the weak boson fusion process or of longitudinal 
weak boson scattering becomes particularly important. 

The study of the "gold-plated" ZZ I'^i^t^t' Higgs boson decay mode alone will 
clearly not suffice to completely understand the nature of electroweak symmetry breaking or 
may not yield a significant signal 0]. Because of severe background problems one needs to 
utilize as much information as possible about the full structure of the underlying qq — >■ qqVV 
process when trying to utilize other channels like W^W~ or hadronic decays of 

the produced weak bosons. This includes forward-jet-tagging to suppress processes like 
qq VV, central-jet- vetoing for the reduction of top-quark backgrounds [Q, looking for 
relatively low overall hadronic multiplicities or searching for rapidity gap signatures 
to exploit the color singlet exchange between the two initial state quarks as a distinguishing 
feature of the signal. In all these cases the knowledge of QCD emission corrections is 
important or even crucial to assess the acceptance of the various procedures for signal events. 

In this paper we present first results of the calculation of all 0{as) real emission correc- 
tions to the full electroweak process qq qqWW (and crossing related ones) with subsequent 
leptonic W decays. Results are obtained by the numerical evaluation of polarization ampli- 
tudes and were implemented for pp and pp collisions. An outline of the calculation is 
given in Section II, details are relegated to an Appendix. 

In two subsequent Sections we apply this calculation to two of the issues mentioned 
above. The suppression of QCD radiation into the central rapidity region is investigated in 
Section HI. Color coherence between initial and final state radiation off each of the incoming 
quarks and the absence of color transfer between the two quark lines in the qq — >■ qqWW 
process leads to strongly suppressed radiation into the rapidity range between the two final 
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state quarks. This confirms the expectation that rapidity gaps, i.e. regions of very low 
hadronic activity, may form between the two quark jets. The influence of processes with 
gluons in the initial state {e.g. gq —>■ qqqWW) on the radiation pattern and the effect of 
radiation on the average gap width are analyzed here. 

In Section IV we then investigate the effects of QCD radiation on forward-jet-tagging and 
central-jet- vetoing as discussed in Ref. [|] . Because of the additional jet activity the efficiency 
of a central-jet- veto might be severely affected by higher order QCD corrections. We show 
that this is indeed the case for the electroweak background of transverse W production, 
but not for the heavy Higgs signal itself. We calculate the signal acceptance of the tagging 
and vetoing techniques and find, typically, a 15% reduction of signal rates compared to the 
lowest order results. 



II. CALCULATIONAL METHODS 

The processes which need to be considered are the QCD real emission corrections to weak 
boson scattering. In the on-shell approximation for the produced W^s this corresponds to 
the cross section evaluation at order 0{as(yQED)j ^^^^ level, of the complete process 

gig2 ^ w^w- g (1) 

and all crossing related processes, like, e.g., 

gqi q^mi w^w~ . (2) 

Some representative Feynman graphs are shown in Fig. |l] and the full set of contributing 
graphs is outlined in the Appendix. We are primarily interested in the weak boson scattering 
process as depicted in Figs, la and b. However, when trying to calculate the distributions 
of the final state quarks we need to treat the incoming weak bosons as off-shell particles. 
Electromagnetic gauge invariance then requires to consider the t-channel photon exchange of 
Fig. lb together with W bremsstrahlung off the two quark lines. We are thus lead to consider 



3 



the full 0{asaQ^jj) process, including all W bremsstrahlung processes off the external quark 
lines, as shown, for example, in Fig. Ic. 

Questions of electromagnetic gauge invariance do not arise when considering the s- 
channel Higgs resonance only, as depicted in Fig. la. For a simplified definition of the 
Higgs boson signal we shall also consider these "s-channel resonance" contributions in iso- 
lation. It is well known, however, that for a large Higgs boson width, i.e. in the case of a 
heavy Higgs boson, the resonance contribution overestimates the production cross section 
of longitudinal W^s at large values of the W pair invariant mass [l^. The "s-channel reso- 
nance" contributions should therefore only be taken as a qualitative estimate of the Higgs 
boson signal. As an alternative means of isolating the effects of a heavy Higgs boson or of 
the scattering of longitudinal weak bosons we shall also consider the cross section difference 
<^{f^H) ~ <^{f^H = 100 GeV), i.e. the excess events over a light Higgs boson scenario. 

The lowest order process, without the gluon emission considered here, has been calculated 
by several groups p!0|J^. We obtain the real emission QCD corrections by numerically 
evaluating the amplitudes, using the helicity amplitude calculus of Ref. as implemented 
in the HELAS package ||^. Fermion masses are neglected everywhere. When considering 
leptonic I^-decays the W propagator factors are taken in the narrow width approximation 
in order to eliminate additional Feynman graphs which are required by gauge invariance for 
off-shell W^s. Details are given in the Appendix. 

A simple Breit-Wigner resonance form with a constant width in the timelike region is 
used for the Higgs boson propagator. Refinements can be made [|lT|, but they are not 



relevant in the following since we are primarily interested here in the QCD structure of 
weak boson fusion processes and do not aim at a precise modeling of the SM prediction for 
W^W~ + 3 jet production. Finally, the phase space integrals over the squared amplitudes 



are performed with the VEGAS integration package |12]. For the parton densities inside 



the proton we use set D'_ of Ref. |]T3[ and we choose the geometric mean of the final state 
parton transverse momenta as the scale of the structure functions and of the strong coupling 
constant asiQ"^). 



Because we are working in the Born approximation and since we are using massless 
quarks throughout, the total cross section for pp — > W^W~X comes out to be divergent 
and we need to impose acceptance cuts to obtain finite results. 

t-channel photon exchange as shown in Fig. Ic is singular at low and in addition 
the parton model no longer provides an adequate description in this phase space region. 
Rather, one should use measured electromagnetic form factors of the proton. While a proper 
treatment of low photon exchange is possible [jl4| one finds that this region can safely be 



neglected when considering signatures for weak boson scattering |^ . Hence, we only consider 
the deep inelastic scattering region in the following by imposing the cut > 4 GeV^. 

Annihilation diagrams contribute in processes with identical fiavors on the two quark 
lines. An example is shown in Fig. Id. Here the splitting of the s-channel photon into a 
qq pair diverges at small dijet invariant masses. Again this process is unimportant for the 
weak boson scattering regime in which we are interested. We eliminate the singularity by 
requiring irijj > 10 GeV for all pairs of final state partons. Finally, collinear singularities, 
which arise from initial state gluon radiation, are avoided by considering partons of finite 
transverse momentum only. 

These cuts eliminate all numerical divergencies. They are largely superseded by the ac- 
ceptance cuts which will be made in actual experiments. Throughout the following we 
require that all pairs of partons which satisfy the jet-identification requirements (typi- 
cally pt > 40 GeV and a limited range in pseudorapidity) are well separated in the 
pseudorapidity-azimuthal angle plane, 

% = (Ar7|^. + A02^.)^ >0.7, (3) 

and similarly that all decay leptons are isolated, 

R,, = {Ar^l + A^iy^ > 0.7 . (4) 

Additional requirements on jet and lepton transverse momenta and on the angular accep- 
tance will be listed separately in the following sections. 



We will want to study the effects of the 0{as) QCD corrections on one and two jet inclu- 
sive distributions and we would like to use the code to assess the probability of soft radiation 
into the rapidity interval defined by two tagging jets. A complete calculation of such observ- 
ables would require the determination of virtual effects and the resummation of soft gluon 
emission, tasks which are clearly beyond the scope of the present paper. Instead we use the 
"truncated shower approximation" (TSA) to obtain estimates of these observables . 

In the TSA we replace the tree level differential cross section for three parton final states, 
da{WWjjj)Tu by 

da{WWjjj)TSA = da{WWjjj)TL (l - e-^'-'-^-Z^?-) , (5) 

Here PTj,min is the smallest transverse momentum of the three final state partons. As 
PTj the final factor in Eq. acts as a regulator of the small pt singularity. The TSA 
parameter ptsa is chosen to reproduce the lowest order cross section within a given set of 
acceptance cuts. For the full electroweak process the singularities associated with low 
photon exchange would introduce a very strong dependence of ptsa on the phase space region 
which is considered. This is not the case, however, for the Higgs boson signal as defined by 
the "s-channel resonance" contribution. It is well known that the i^-factor for the qq qqH 



process is close to unity llg] and this process has a finite total cross section at tree level. By 
choosing Ptsa such that the "s-channel resonance" contribution to (y{WW jij)^^,^ reproduces 
the corresponding lowest order cross section, we obtain an algorithm which agrees with the 
full 0{a 

sO^qed) calculation in the 3 jet phase space region, which allows to study 1-jet and 2- 
jet inclusive distributions to the same order, and which provides an excellent approximation 
to these distributions and to the overall normalization at full Olasag^j;)) in the phase space 
region where the Higgs resonance dominates. 

Anticipating the acceptance cuts to be used below, we have determined the TSA param- 
eter by matching the total Higgs production cross sections in pp scattering at ^/s = 40 TeV 
with one or two visible jets exceeding a minimal transverse momentum. Results are given in 
Table I. One finds that ptsa is quite insensitive to the value of the Higgs boson mass and also 
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varies very little with the pTj,min cut imposed in the 1-jet inclusive case. When requiring 
the presence of two partons of px > 20, 40, or 60 GeV this variation is somewhat stronger, 
but even here the variation in ptsa of ~ ±1 GeV corresponds to normalization changes of 
the cross section by up to 30% only, which is within the uncertainty range of our tree level 
calculation. 

III. COLOR STRUCTURE AND RAPIDITY GAPS 

The characteristic features of the 0{as) emission corrections to weak boson scattering can 
be understood in terms of the color structure of the lowest order process. In the dominant 
t-channel contributions (see Figs, la-c) no color is exchanged between the two incoming 
quarks (or anti-quarks). For small scattering angles of these incoming fermions the color 
charges are accelerated rather little and the resulting "synchrotron radiation" of gluons 
occurs predominantly into the forward direction, between the beam axis and the direction 
of the scattered (anti) quark. One thus expects very little hadronic activity to arise from 
bremsstrahlung off the hard process in the central region, between the two scattered quarks. 
Since the energy carried by the quarks is typically much larger than the virtuality of the 
emitted weak bosons (of 0{mw)) these quarks are produced at rather large pseudorapidities, 
leaving a wide rapidity region with suppressed radiation 0. This pattern is fundamentally 
different from the one expected for typical background processes (like e.g. gg —>■ tt ^ 
bbW~^W~) where color exchange between the two incoming partons leads to strong gluon 
radiation into all regions of the legoplot. 

The suppressed radiation into the central region may then lead to the formation of 
rapidity gaps, pseudorapidity regions of very low or no hadronic activity except for the 
Higgs decay products This rapidity gap may be filled, however, by hadrons from the 

underlying event ||T^ and the full formation of the gap is expected to be observable only 
in a small fraction of all signal events where fluctuations or the absence of multiple parton 
interactions make the radiation pattern of the hard scattering process visible []^,|TB|,|TU[ . 
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Here we do not discuss this last question of the gap survival probability any further, rather 
we investigate whether the higher order QCD corrections indeed do lead to the radiation 
patterns suggested by the color flow arguments. Processes with three quarks in the final 
state (like .e.g gq qqqH) might substantially alter the expected pattern. In addition, a 
precise knowledge of the radiation in signal vs. background events may help to distinguish 
the two even in the presence of an underlying event [^]. 

The features of the heavy Higgs boson signal discussed above suggest a search following 



the strategy developed for rapidity gap events in Ref. pOf. We study 2 jet inclusive WW 
events where both W^s decay leptonically and where the two charged decay leptons fall into 
the central rapidity region between the two tagging jets. More precisely we start with events 
containing two charged leptons of high transverse momentum in the central region, 

PTi > 100 GeV, \r],\ < 2 . (6) 

No hadronic jet with transverse momentum 

PTj > 40 GeV (7) 

is allowed in the pseudorapidity interval between the two leptons. On either side of the 
lepton pair we search for the first jet of ptj > 40 GeV, i.e. we require the presence of two 
tagging jets. The pseudorapidity range between the tangents to the jet definition cones (of 
radius R = 0.7) is called the "gap region" and the two charged leptons are required to fall 
into this gap region. Denoting by rjj-^ < rjj^ the pseudorapidities of the two tagging jets we 
thus require 

r]j^ + 0.7 < rji- < rjj^ - 0.7 . (8) 

The width of the gap region is denoted by 

Vgap = \Vji - Vj2\ - 2 ■ 0.7 . (9) 

The average gap width is expected to be large for weak boson scattering events but not 
necessarily in the phase space region where W bremsstrahlung dominates. This is confirmed 
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in Fig. where da/di/gap, within the acceptance requirements of Eqs. ®-(|]), is shown 
for the production of a rriH = 800 GeV Higgs boson in pp colhsions at y/s = 40 TeV. 
The Higgs boson and longitudinal weak boson scattering contributions can be isolated by 
comparing with the case of a very light Higgs boson, taken here as mn = 100 GeV. For 
the lepton acceptance requirements of Eq. (|^) the resonance contribution of the light Higgs 
boson is completely negligible. Thus the rriH = 100 GeV curve represents an estimate of the 
contributions due to W bremsstrahlung and transverse W production, which may be called 
an "electroweak background" to the heavy Higgs boson signal. Comparison of the light and 
heavy Higgs boson curves in Fig. ^indeed shows that the signal events are produced at large 
average i/gap. 

The main effect of the third parton in the final state is a slight narrowing of the average 
gap width from < i/gap >= 4.9 in the lowest order calculation to < i/gap >= 4.5 in the full 
0{as) result (for rriH = 800 GeV). There are at least two reasons for this effect: the possible 
emission of a third jet into the region between the original two quark jets redefines the gap 
region and in addition the emission of a third parton outside the gap region, i.e. at very large 
rapidities, raises the center of mass energy of the event and hence leads to a kinematical 
reduction of the gap width. Neither effect is very large, however. 

The QCD radiation of the third parton is dominated by collinear emission close to one 
of the two tagging jets. Because of the coherence between initial and final state radiation 
the emission occurs mainly between the tagging jets and the beam axes, i.e. outside the gap 
region. This effect is demonstrated in Fig. |^ where the distribution in the rapidity difference 
of the third (typically soft) parton and the closest tagging jet is shown, 

Arjsj = sign ■ (?7soft - VjAosest) ■ (10) 

The sign is chosen such that positive Arjgj corresponds to radiation outside the interval 
marked by the two tagging jets and negative values correspond to radiation towards the gap 
region. The edge of the gap region is visible as a step at Arj^j = —0.7 which is a result of 
the jet separation requirement of Eq. (^. Clearly emission of the third parton outside the 
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gap region is preferred. This is true both for the heavy Higgs case (sohd hne) and for the 
electroweak transverse W background, simulated by the uih = 100 GeV scenario (dashed 
hne) . 

The resulting radiation pattern is best appreciated by choosing the center of the gap 
region as the origin and then normalizing the pseudorapidity of the third parton to the gap 
width. This is achieved by using the variable 

l^ji -^j2l -1-4 ■ 

Thus z = corresponds to to the gap center while z = ±1 indicates the edges of the gap 
region. The resulting radiation pattern, da/dz, is shown in Fig. The probability for 
radiation into the gap region is strongly suppressed both for the heavy Higgs signal and for 
the electroweak background. The radiation pattern which is expected for t-channel color 
singlet exchange is thus confirmed at 0{as). 



IV. QCD EFFECTS ON JET-TAGGING AND CENTRAL-JET- VETOING 

When trying to observe the W^W~ decay mode of a heavy Higgs boson at a hadron 
collider, one needs to fight serious physics backgrounds. The production of a top-quark pair 
with subsequent decay ti hW^hW~ is particularly troublesome due to the large rate of 
this background. Of somewhat lesser importance is qq —>■ W^W~ pair production ||^. Higgs 
boson production via weak boson fusion differs in important aspects from these background 
processes. Forward-jet-tagging of the quark-jet (s) off which the initial state W^s or Z's were 
radiated offers a powerful tool for background suppression [0. In addition one can use the 
fact that the VT's arising from top-quark decay are always accompanied by a nearby 6-quark 
which often manifests itself as an additional hadronic jet. When searching for the leptonic 
decays of the produced W^s in the central rapidity region a veto on any additional central 
jet constitutes a powerful tool for top-background reduction. 

Single forward-jet-tagging and central-jet- vetoing were analysed in Ref. |Q. The Higgs 
boson signal was simulated with a full 0(aQED) Monte Carlo program. While this appears 
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adequate for the analysis of tagging jet distributions, the acceptance of signal events under 
severe central-jet- vetoing conditions might be strongly degraded by 0{as) QCD corrections. 
This question can be addressed with the tools described in Section 11. 

Following Ref. ^ we consider purely leptonic decays of the two W bosons and concentrate 
on the Higgs resonance region by imposing stringent lepton acceptance cuts, 

\r]i\ < 2 , pTi> 100 GeV, Apru > 400GeV . (12) 

Here ApTu = |pt£i — Pt£2\ is the difference of the two charged lepton transverse mo- 
menta [p^] . 

In a second step one requires the presence of at least one hadronic jet of transverse 
momentum 

PTj > 40 GeV . (13) 

The most energetic such jet is the tagging jet candidate. In Figs. ^ and || the pseudorapidity 
and transverse momentum distributions of this jet, as obtained with the full C^(«s«qed) 
Monte Carlo, are shown for a. mu = 800 GeV signal and the electroweak background as 
defined by a mu = 100 GeV scenario. Most of the tagging jet candidates at low r7j,tag arise 
from transverse W bremsstrahlung and hence a stringent rapidity cut on the tagging jet is 
highly efficient for the heavy Higgs signal. The transverse momentum of the tagging jet in 
Higgs signal events is typically quite low, while the electroweak background distribution is 
rather fiat, reflecting the dominance of W bremsstrahlung: the strong lepton p^ cut selects 
high pt VT-bosons which were radiated off high transverse momentum quarks. The tagging 
jet distributions are quite similar to the 0{aq^j^) calculation which, therefore, are not shown 
separately. As in Ref. 0] signal events are selected by requiring the presence of a tagging 
jet of 

E,,tag > 1 TeV , 3 < |r/,,tag| < 5 . (14) 

Because of the presence of a third colored parton in the final state the full 0{as(y.q^j)) 
calculation yields a larger probability to find an additional jet in the central region than 
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the corresponding result without QCD corrections. Let us define the veto jet candidate as 
the remaining parton, other than the tagging jet, with largest transverse momentum. The 
pseudorapidity distribution of this veto jet candidate is shown in Fig. |^. 
In Ref. m a veto of any jet in the central region satisfying the conditions 

\Vjiveto)\ < 3 , pT(veto) > 30 GeV (15) 

was found to provide adequate top background rejection while retaining most of the signal 
events when neglecting QCD corrections to the Higgs signal. The cross sections for the 
Higgs boson signal after vetoing central jets with varying minimum transverse momentum 
requirements are shown in Table II for both the lowest order "2 parton calculation" and 
the 0{as) "3 parton calculation". Also given are the efficiencies of retaining the Higgs 
signal after central-jet-vetoing. The 100% level corresponds to the cross section difference 
aijUH = 800 GeV) — (j{mH = 100 GeV) before rejecting events with jets of |?7j| < 3 and 
transverse momenta as listed in the table. 

In the phase space region which is dominated by W bremsstrahlung (as simulated by 
the ruH = 100 GeV scenario) QCD radiation does indeed strongly affect the efficiency of 
a central-jet-veto. In fact this leads to a further reduction of the electroweak background 
compared to the Higgs boson signal. At 0{as) the signal acceptance of the central-jet- veto 
is lowered by ^ 20% and the resulting signal cross sections are typically only ^ 15% lower 
than before QCD corrections. Hence the conclusions of Ref. remain valid when taking 
real emission QCD corrections into account. 
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APPENDIX A: MATRIX ELEMENTS 



The external state of the processes considered in Section |I| consists of four quarks, one 
gluon, a and a W~ boson. The bosons will be taken later to decay into leptons, 
but this does not change the nature of the calculation. The four (anti) quarks are labelled 
as a, b, c, d, along with their corresponding momenta pa,Pb,Pc,Pd, and helicities Xa, Xb, Ac, A^. 
The four (anti) quarks are paired to form two fermion lines, denoted as (ac) for fermion 
flow from a — >■ c, and [bd) for fermion flow from b ^ d. The gluon is taken to have a 
four- momentum given by Pe, and helicity Ag. The momenta of the and W~ final-state 
bosons are denoted by qi and q2, respectively. 

The qqqqW^W~ g amplitudes are evaluated numerically using the HELAS program pack- 
age 0. The HELAS program contains subroutines evaluating and multiplying the various 
factors in the Feynman graphs. Hence we shall employ a notation directly relating to the 
structure of HELAS. To begin, the generalized propagator for a vector boson of mass M in 
the unitary gauge is given by 



D>"'{q^, M) ^ 



^ M2 



(Al) 



g2-M2 + ^Mre(g2) 

with q^ representing the momentum transfer squared, and F denoting the q^ independent 
decay width of the intermediate particle. The step function B(g2) eliminates the imaginary 
part of the vector boson propagator for spacelike momentum transfer. The analogous choice 
for the heavy Higgs boson propagator will introduce only a small error compared to unitarity 
requirements Similarly, for massless particles, the propagator function is defined as 



d'^iq^) = --^r (A2) 

r 

Vertex insertions are defined via 

yt = r [g'v + 9\7'] (A3) 

where the index i on the vector and axial coupling constants gy and g\ will list the vector 
bosons and fermions coupling at a particular vertex. 
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The zero- width- approximation is employed in order to include spin-correlation effects of 
the decay products. The polarization vectors of the final state W bosons are denoted 
by el{qi) for the boson and £2(^2) for the W~ boson. 



Ai)* = JUki,k2,PuP2)JjT^ (A4) 

V ^vIwL W 



4{Q2, A2)* = -^31(^3, h, P3, Pa)\ -irrir- (A5) 

V JyJ-wi- w 

where the truncated decay current is defined in terms of the final state fermion spinors 

Jf2 = u{ki, pi)V-^ffv{k2, P2) (A6) 
^34 = ^^(^3, Pz)V^ffv{ki, P4) (A7) 

Inserting these decay currents for the W polarization vectors the resulting matrix elements 
will describe the full process 

QaQc ^ qbQdW^W-g , W+ ^ 44 , W- ^ i^h , (A8) 

and all processes related by crossing of the quarks and gluons. 

In order to simplify the spinor algebra we use a bra and ket notation: 





(A9) 


u{pb,Xb) = \b) 


(AlO) 


HPc, Ac) = (c| 


(All) 


u{Pd, Ad) = {d\ 


(A12) 



and similarly, the wavefunctions for antiquarks are defined as 



v{Pa: -Aa) = {a\ 


(A13) 


v{Pb, -Aft) = {b\ 


(A14) 


v{Pc, -Ac) = |c) 


(A15) 


v{pd, -Ad) = \d) 


(A16) 
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The currents, corresponding to the sphtting of a virtual 7, Z into a final-state W^W 
pair via the triple-boson vertex, are 



(A17) 
(A18) 



where E represents the WW'f triple-boson vertex rule, with the index ordering corresponding 
to incoming W^W~j bosons or outgoing W~W~^j bosons. The tensor F represents the 
WWZ triple-boson vertex rule in analogous notation. 

Off^-shell fcrmion wavef unctions, corresponding to gauge boson emission off the fermion 
lines, are defined via 



|l,a) 


Pa- i\ 


(A19) 


|2,l,a) 




(A20) 


\U,a) 




(A21) 


(c,2| 




(A22) 


(c,2,l| 




(A23) 


{c,U\ 




(A24) 



Similar expressions hold for wavef unctions involving fermions h and d. 
Gluon insertions along a fermion line are notated as 



\g, l,a 
1^,2, 1, a 



AX* 



i 



i 

I 



"a 



i>a- i>e- i\- ii 



efV;^^^|2,l,a) 



i 



{c,g\ = {c\V^^fe^^ 
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(A25) 
(A26) 
(A27) 
(A28) 
(A29) 
(A30) 



c,^,2| = {c,g\V^fUT ^ , ' , , (A31) 



i 



(c.2,,,l|Hc.2,9|V;"'".r^-^;^ (A33) 

(c.2.1,.|^(c.2.1|vr.r ,^^^;^^^ (A34) 
Simple electroweak currents are denoted as in 

Jt;{ac) = #"((Pa -Pc)')(c|K/^|a) (A35) 

J{^(ac) = ^'^''((p^ -p,)^M^)(c|K'^^^|a) (A36) 

J^(ac) = D^'^iip, - p,f, Mz){cK'"\a) (A37) 

J^{ac) = cos J^(ac) + sin OwJUjiac) (A38) 

with (ac) representing fermion current flow from a ^ c. Fermion current flow from c — > a 
would be represented by the argument (ca). 
When including gluon emission the notation 

J^^iac) = d^-iiPa-Pc-Pe?){c\Vjff\g,a) (A39) 

Ji;^{ac) = d^'^iipa -Pc- Pe)'){c, g\Vjff\a) (A40) 

J^^ica) = d^''{{Pc-Pa-Pef){a,g\Vjf^c) (A41) 

Ji}^(ca) = d^''{{p,-p,-pe)'){a\V;j"\g,c) (A42) 

is used, and similarly more complicated single W boson emission currents are denoted by 

J^^iac) ^d^''{{p,-p,-q^f){c\Vj"\l,a) (A43) 

J^^iac) = D^mPa-pc-qi)\Mw){cK'^^^\l,a) (A44) 

J^,{ac) ^ D^''{{pa-pc-qi)',Mz){c,lK^"\a) . (A45) 

Including the emission of a gluon, the notation of the latter three changes to 

J,V(ac) = d^'^iipa-Pc-Pe - qi)'){c, l\Vj"\g,a) (A46) 

16 



Ji^^,{ac) = D^^iipa -Pc-Pe- qi)\ Mw){c, 1, ^I^T^^Ia) (A47) 
JzigM = D^''i{pa-p,-p, - q,)\Mz){c,g, lK^"\a) (A48) 

and similarly for interchange of outgoing bosons 1 and 2, interchange of fermion current flow 
a c, interchange of fermion line (ac) for {hd), or a combination of any and/or all of the 
above interchanges. Notice that in all cases the space-time index of the current corresponds 
to the explicitly listed electroweak boson in the index identifying the current. 

Double W boson emission currents may be defined similarly, as in the examples 



J^,^{ac) = d'^^'iiPa-Pc -qi- q2f){c\Vj"\2, 1, a) (A49) 

J(V2(«c) = Z?^"((pa -Pc-qi- q2)\ Mw){c, 2|K^^^|1, a) (A50) 

J^^(ac) = D^^iipa -Pc-qi- q2)\ Mz){c, U\V,^"\a) . (A51) 

Including gluon emission another index needs to be added, 

J^,,^{ac) - d^'^ipa -Pc-Pe-qi- q2f){c, 2\Vj"\g, 1, a) (A52) 

J^guiac) = D^'^iipa -p^-p^-q,- q^f ^ Mw){c, U,g\V;^"\a) (A53) 

Jgzuiac) - D'^^iipa -Pc-Pe-qi- ^2)', Mz){c, 2, lK^"\g, a) , (A54) 



where again, all possible interchanges of outgoing bosons, fermion current flows, and fermion 
lines are allowed. 

For sums of currents involving gluon emission off the same fermion line, a semi-colon 
notation is introduced. For example, 

4,(ac) - J^^iac) + Jl^giac) (A55) 

Jzi2;giac) = J'^zA^c) + JzgAo^c) + ■^Zlg2(ac) + Jzi2g{o^c) (^56) 

represent a photon and a Z current given by the sum of currents involving a gluon emission 
along all possible topologically distinct locations on the fermion line. In terms of the bra 
and ket notation, the semi-colon notation is analogously introduced by example. 
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(c, 1, 2|X|a).^ = (c, 1, 2\X\g, a) + (c, 1, 2,g\X\a) 

+ (c,l,(7,2|X|a) + (c,(7,l,2|X|a), (A57) 

where X is an arbitrary vertex insertion. 

The scalar current obtained from WW fusion into a Higgs boson is defined via 

HWWh ^ ^((^^ ^ Mh)/,^^^ (A58) 

where I^^^ represents the WWH vertex rule. Similarly, 1^,^^ represents the ZZH vertex 
rule. 

The SU{2) quadruple-boson vertices are represented by S and T. S^i^xp represents the 
WWWW effective coupling, which includes the contact term, as well as s-channel and t- 
channel 7, Z exchange terms. T^j^Ap represents the WW^WW^ effective coupling, which 
includes the contact term, as well as t-channel and w-channel W^ exchange terms. 

The tensor G represents the WWW^ triple-boson vertex rule, with the index ordering 
corresponding to incoming W~^W~W^ bosons, or outgoing W~W~^W^ bosons, respectively. 
The current J3 is, as defined in Eq. (ASS), a SU{2) superposition of 7, Z currents. 

With this notation we can now express the amplitudes, which correspond to the various 
Feynman graphs, in a very compact notation which can directly be translated into a calling 
sequence for HELAS subroutines. For the Feynman graphs of Figure |.a one gets 

- ^M,ia) = Ji^Jac)J^ibd)lZ'^^HZ'^''e',*eT (A59) 



~tM,ib) = J^JM)J(;.(ac)/™i/,^^M*er (A60) 



where Ai{a) represents an amplitude including all distinguishable gluon emission insertions 
along the fermion line a — > c. Similarly, A4{b) represents an amplitude including all distin- 
guishable gluon emission insertions along the fermion line b ^ d. 

Only the Ai{a) amplitudes are shown below, with the A4{b) amplitudes derived via 
interchange of the fermion indices (ac —>■ bd). The remaining M.[a) amplitudes can be 
represented as follows 



up 



-1 "=2 



(A61) 
(A62) 



with corresponding Feynman diagrams presented in Figure ||.b. Notice that for a given set 
of external quark flavors only one of the two amplitudes Ai2 or will be nonvanishing: 
A^2 corresponds to a coupling to the (ac) quark line while for Ai^ the emission of a 
W~ off the (ac) quark line must be allowed. Analogous flavor selection rules will be implied 
for all remaining diagrams as well. The amplitudes for these are given by 



— iM4^{a) 






(A63) 








(A64) 


-iMe{a) 




J^^,(ac)J^(M)/,^/^if™e^er 


(A65) 


—iMjia) 




J3^,(ac)J3^(M)T,,,.6^6r 


(A66) 


onding Feynman diagrams presented in Figure 9 




-iMs{a) 




J^W;9MJ^ibd)G,xA* 


(A67) 


—iMgi^a) 




J^2,g{cic)r^ipd)G,xA* 


(A68) 


—iMio{a) 




J^w,g{(^c)r,ipd)Gx,A* 


(A69) 


-iMn{a) 




Jwi;9iac)J^{bd)Gx,ue'2* 


(A70) 


-iMuia) 




J^.g{ac) {bd) Gyx^ie^* 


(ATI) 


-iMisia) 






(A72) 


— iA^i4(a) 




Jl^.g (ac) JiV {bd)Gxu^le2* 


(A73) 


-iA^i5(a) 




J^.g (ac) J(^i {bd)Gxu^ie2* 


(A74) 


—iMie{a) 




Jw;g{ac)r2z{bd)G^xA* 


(A75) 


-iMiria) 




Jw;9{ac)JUhd)G^Xue'r 


(A76) 


-iMisia) 




Jw;9{ac)JUhd)Gx^,A* 


(A77) 


—iA^ 19(a) 




Jw;gi(^c)J!^,{bd)Gx,A* 


(A78) 


-iM2o{a) 




J^3;gMJ!^{bd)G,xA* 


(A79) 



19 



~iM2i{a) 


= J^2;gi^c)J!^{bd)G^x^,e^* 


(A80) 


-iM22{a) 




(A81) 


-iM23ia) 


= 73^1.3 (ac)J(;,(M)G'A.^e^* 


(A82) 


— ■iA^24(a.) 


= J^2rA^c){d\V^^f\b) + J^2Z-,g{ac){d\V^^^b) 


(ASS) 


-iAl25(a) 


tH, / \ / 7 1 T 7-'*v f f \ 1 \ tII, / \ 1 i\Tr^ff\i\ 

= J^^2;g{^c){d\V;:^f\b) + Jfz2;g(ac)(rf|\//^^|6) 


(A84) 


— iA^26(a) 


= JiJi2;g{ac){d\V^"\b) + J^i2-giac){d\V;^^b) 


(A85) 


—iM27ia) 


= -^2l7;9(«c)(rf|^7^|fe) + J21Z;giac){d\V;^^b) 


(A86) 


-ZA^ 28(a) 


= J2yl;g{''c){d\V;:^f\b) + J^zi-,giac){d\V;^^\b) 


(A87) 


-iTW 29(a) 


= ^721;^?(«c)(rf|\/7^|6) + J^21;gM{d\V;^f\b) 


(A88) 


-zA^ 30(a) 


= J'i2^{bd){c\V;ff\a).g + Jr2z(M)(c|\//^^|a).^ 


(A89) 


— zA^3i(a) 


t// /7 7\ / It 7''V f f \ \ . tII /I 7\ / \ T r 7, f f \ \ 

= J{;2(M)(c|\//^|a).^ + J^z2{bd){c\V^ff\a).g 


(A90) 


-zA^ 32(a) 


~rll / 7 7\ / It 7"V f f \ \ , ~rll /7 7\/ iTrJ^ffl \ 

= Ji^i2ibd){c\V;jff\a).g + J^i2ibd){c\V;"\a).g 


(A91) 


-?Al33(a) 


= J2i^ibd){c\V;^^f\a).g + J^izibd){c\Vff^\a).g 


(A92) 


—zAl 34(a) 


= J2V(^^)(c|^/^l«);, + J^ziibd){c\V^'"\a).^^ 


(A93) 


-iX35(a) 


= J,^2i(M)(c|^/^|a)^^ + rz2iibd){c\V^'"\a).^^ 


(A94) 


onding Feynman diagrams presented in Figure IC 




-zAl 36(a) 


= Jr,^,(ac)(d|y//|2,6) + J^,.^^{ac){d\V^^^f\2,b) 


(A95) 


—iAI 37(a) 


= J!;,.^g{ac){d\v;!ff\2,b) + rz,.^^{acmv^ff%b) 


(A96) 


-zAl 38(a) 


= J,%g{ac){d,2\V;j^f\b) + J^,.^^{ac){d,2\V^ff\b) 


(A97) 


-iAl39(a) 


= Ji;,,g{ac){d,2\V;j^f\b) + r,,.^^{ac)MV^ff\b) 


(A98) 


— zAl4o(a) 


= J,%g{ac){d\V;^^f\l,b) + J^,.^^{ac){d\V^ff\l,b) 


(A99) 


— zAl4i(a) 


= Ji;,,g{ac){d\v;^^f\i,b) + r,,.^^{ac){d\v^ff\i,b) 


(AlOO) 


—zAl 42(a) 


= J2%,{ac){d, l\V;j^^b) + J^z;g{<^c){d, l\V^'^^b) 


(AlOl) 


— iAl43(a) 


= J,\,(ac)(rf, 1|V7^^|6) + rz2;g{ac){d, 


(A102) 


—zAl 44(a) 


= Ji}r,,{ac){d\V^ff\b) + Jf^z;g{ac){d\V'ff\b) 


(A103) 
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a) 




(A104) 




a) 


= 'Ju^ibd){c\V^f^a).^ + Ji^zipd){c\V^ff\a).g 


(A105) 




a) 


= J!;jj{hd){c\V^'ff\a).^ + J''zuipd){c\V^ff\a).g 


(A106) 


-iMis 


a) 




(A107) 


-iMaq 


a) 


tU / \ / 7 1 T 7-14/ f f \ 1 \ 

= J21W;giac){d\V^^f\b) 


(A108) 




a) 


tU / \ / 7 1 T f f \ 1 \ 

= Jw21;giac){d\V^^f\b) 


(A109) 




a) 


= Jwi2;giac){d\V^^^b) 


(Alio) 


-iM52 


a) 


tU /I 7\ / \ T rW f f \ \ 

= J2iwibd){c\V^^^\a).g 


(AUl) 




a) 


= •Ji2wibd){c\V/;'^^\a).g 


(A112) 




a) 




(A113) 


-iM55 


a) 




(AIM) 


ending 


Feynman diagrams presented in Figure 11, and finally 






a) 


= J^W;giac){d\Vi^^f\l,b) 


(A115) 


-iM57 


a) 


= J{^2;giac){d\V^ff\l,b) 


(A116) 




a) 


= J^W;giac){dA\V;^f^\b) 


(A117) 


-iMm 


a) 


= J{^2;giac){dMV;^ff\b) 


(A118) 


-iMqo 


a) 


= J^W;giac){d\V^f^\2,b) 


(A119) 




a) 


= J{^,,,{ac){d\V^f^2,b) 


(A120) 


-iMm 


a) 


= ^iV;.(«c)(rf,2|v;^//|6) 


(A121) 


-iMe3 


a) 


= J(^.,^,(ac)(d,2|V;^//|6) 


(A122) 


-iMqa 


a) 


= Ji}w;gM{d\V^ff\b) 


(A123) 




a) 


= J^U;giac){d\V^ff\b) 


(A124) 


-iMm 


a) 


= Ji;wibd){c\V^"\a),, 


(A125) 




a) 


= J{^uibd){c\V^"\a).^^ 


(A126) 



with corresponding Feynman diagrams presented in Figure |I2|. Notice that there are no 
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diagrams involving Goldstone bosons. This is a remnant of the choice of the unitary gauge 
for this calculation. 

Once the individual amplitudes are evaluated, the two amplitude sums are calculated as 

Mia) = Miia) + Miia) + ...+ Meria) (A127) 
M{b) = Miib) + M2ib) + ... + Merib) (A128) 

The overall amplitude including the QCD color structure is hence 

M = (a) + 6acTi,M (6) (A129) 

where a, b, c and d represent the color indices of the four quarks, and i represents the color 
index of the attached gluon. The required sum of the amplitude square, |A^P, over all 
external colors thus yields 

J:\M\' = iV^^ ( \M{a)\' + \M{b)\' ) (A130) 

colors 

for an SU{N) gauge theory {N = 3 for QCD). 

When identical final-state fermions are present one must anti-symmetrize the result- 
ing matrix element. This complicates the resulting color structure. In comparison with 
Eqn. |A129| above, the resulting amplitude contains both t and w-channel processes as 



M = nMMia) + 5acTLM{b) 

- n^6kcM (a) [c^d]- dadTlM (b) [c ^d] (A131) 



This leads to modified terms in the analogue of Eq. ( |A130| ). Defining 



M[ = M{a)+M{b)-M (a) [c^d]-M (b) [c ^ d] (A132) 

M2 = M{a)-M{b) (A133) 

M'3 = M (a) [c^d]-M (b) [c ^ d] (A134) 

M'^ = M{a) +M{b) +M{a)[c ^ d] +M{b)[c ^ d] (A135) 



one obtains 
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E \M\' = J:c;\M'^\' (A136) 



with 



colors j=l 



C[ = {N + 2) [-^j (A137) 
/N^ -1\ 

C'2^N ^— (A138) 



2 

AT (^^] (A139) 



/ iV2 - 1 \ 

= (TV - 2) J • (A140) 

Within SU (3) one needs to insert the color factors C[ — 20, C2 = Cg = 12 and C4 = 4. Inclu- 
sion of s-channel diagrams, and possible t and s-channel interference proceeds with identical 
color structures as for the i-channel, and t and u-channel interference cases, respectively. 
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TABLES 

TABLE L Values (in GeV) of the ptsa parameter of Eq. (|5|) which are needed to match the cross 
sections at 0{as) with the lowest order results for the "s-channel resonance" process. Values are 
given for two Higgs boson masses and different transverse momentum requirements for "visible" 
jets in pp collisions at ^/s = 40 TeV. 

1 jet inclusive 2 jet inclusive 

TJiH = 500 GeV mn = 800 GeV run = 500 GeV = 800 GeV 

PTj > 20 GeV 5.8 5.7 6.7 6.5 

PTj > 40 GeV 6.1 5.9 7.9 7.6 

PTj > 60 GeV 6.4 6.3 8.9 8.9 
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TABLE II. Cross sections (in fb) and central-jet-veto efficiencies e for pp collisions at 
y/s = 40 TeV. Only leptonic decays of both 1^'s {W^ — > e^v,ijL^v) are considered and the lepton 
and tagging jet cuts of Eqs. (|l^)-(14) are imposed, leaving B ■ (j{mH = 800 GeV) = 8.2(7.6) fb 
and B ■ a{mH = 100 GeV) = 1.9(1.5) fb before central-jet-vetoing for the 0{as) (lowest order) 
calculation. A veto is imposed on central jets of < 3 and transverse momenta pTj > as 
given in the table. The signal is defined as the difference of the mn = 800 GeV and mn = 100 GeV 

cross sections. 

a) 2 parton calculation 
B ■ aimn = 800 GeV) B ■ aimn = WO GeV) S • cjsig„ai e 

4.4 0.5 3.8 0.62 

4.5 0.5 4.0 0.65 
4.7 0.6 4.1 0.67 
5.1 0.6 4.5 0.73 



= 20 GeV 
30 GeV 
40 GeV 
60 GeV 



b) 3 parton calculation 

B ■ a{mH = 800 GeV) B ■ a{mH = WO GeV) S • cjsig„ai e 

3.4 0.2 3.1 0.50 

3.6 0.2 3.4 0.54 

3.8 0.2 3.6 0.57 

4.3 0.3 3.9 0.62 



p5?*° = 20 GeV 



30 GeV 
40 GeV 
60 GeV 
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FIGURES 

FIG. 1. Representative Feynman graphs contributing to the electroweak process 

pp W^W~jjjX. 

FIG. 2. Gap width distribution, da/dygap, in the electroweak process pp — > W~^W^ jj{j)X at 
^/s = 40 TeV. Results are shown for the full 0{asaQ^j-)) calculation and Higgs boson masses of 
niH = 800 GeV (solid line) and uih = 100 GeV (dashed line). The cuts of Eqs. (^)-(^ are imposed, 
including lepton acceptance cuts of pxi > 100 GeV, jry^j < 2. The cross section is normalized to 
include all W decay modes, however. 

FIG. 3. Distribution in Ary^j, the minimum pseudorapidity separation of the third (typically 
soft) parton from the two jets defining the boundary of the gap region. Negative values of Arjgj 
correspond to the emission of a soft parton (px < 40 GeV) into the gap region. Curves are shown 
for the same parameters as in Fig. ^. 

FIG. 4. Normalized gap width distribution, da/dz (see Eq. ([ll|)), for W^W~ jj{j)X events in 
pp collisions at ^/s = 40 TeV. z = corresponds to the center of the gap and z = ±1 mark the gap 
edges. As in Fig. I the cross sections are normalized to include all W decay modes. Results are 
shown for itih = 800 GeV (solid line) and uih = 100 GeV (dashed line) and the cuts of Eqs. (^-(|8|) 
are imposed. 

FIG. 5. Pseudorapidity distribution of the tagging jet candidate (the most energetic jet of 
PTj > 40 GeV) in pp — > W~^W~ j{jj)X events at ^/s = 40 TeV. Curves are shown for 1-jet inclusive 
events at 0{as) and Higgs boson masses of mn = 800 GeV (solid line) and niH = 100 GeV (dashed 



line). The acceptance cuts of Eq. (12) and Eq. ( |13D are imposed and the curves are normalized as 
in Fig. |. 

FIG. 6. Transverse momentum distribution of the tagging jet candidate in 

pp W^W~j{jj)X events at -^s = 40 TeV. Acceptance cuts and normalization of the curves are 
as in Fig. [sj. 
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FIG. 7. Pseudorapidity distribution of the veto-jet candidate in pp — > W^W~j{jj)X events 
at ^/s = 40 TeV. Curves are shown for 1-jet inclusive events at 0{as) and Higgs boson masses of 
niH = 800 GeV (sohd hne) and mn = 100 GeV (dashed hne). In addition to the acceptance cuts 
as imposed in Figs. ^ and ^ the presence of a tagging jet of pTj > 40 GeV, Ej > 1 TeV is required 
in the pseudorapidity range 3 < |?7j,tag| < 5. The curves are normahzed as in Fig. 

FIG. 8. Feynman diagrams corresponding to amphtudes A^i through Ai^. 7W3 is obtained 
from by interchanging the and W~ hues in graph b). The crosses mark the four locations 
where the external gluon must be attached, leading to a total of eight permutations for graph b). 

FIG. 9. Feynman diagrams corresponding to amplitudes A^4 through Mj. The circles repre- 
sent the four-boson- vertices as well as s-, t- and/or li-channel electroweak boson exchange which, 
in HELAS, are combined into a single subroutine. 

FIG. 10. Feynman diagrams corresponding to amplitudes TWg through Ai^^. The permuta- 
tions involve the different attachments of the final state W^s to the lower quark line, the interchange 
of final state and W~ , the interchange of upper and lower quark line, and the different at- 
tachments of the gluon. 

FIG. 11. Feynman diagrams corresponding to amplitudes M^q through Ai^^. In determining 
the number of permutations, 7 and Z exchange are not counted separately. 

FIG. 12. 
Feynman diagrams corresponding to amplitudes Ai^e through Aiej. 
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